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Abstract 

Electronic state modulation of the armchair (6,6) carbon nanotubes filled with a linear assembly of oxygen molecules is addressed 
theoretically. Ferromagnetic coupling of encapsulated oxygen produces a magnetic field with cylindrical symmetry, which deviates 
the electron's eigenenergies from those prior to the oxygen absorption. An intriguing spectral gap arises near the Fermi energy, at 
which the gap formation is allowed only when the tube length equals to a multiple of three in units of carbon hexagon. A possible 
means to detect the selective gap formation is discussed. 
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1. Introduction 

The central cavity of carbon nanotubes is an ideal container 
for atoms and small molecules. Initiated by the massive H2- 
storage observation [1], many attempts to introduce foreign en- 
tities into the hollow space of nanotubes have been reported: 
H 2 S H, H 2 2 ll, C0 2 H, C 2 H 2 01, and ArlH are 
only a few to mention. In addition to the potential utility for 
energy-storage materials and molecular transport devices, the 
filling of carbon nanotubes is believed to open up novel appli- 
cations of drug delivery to the ce ll 0, [T^]. In light to me- 
chanical flexibility of nanotubes [11], furthermore, effects of 
their cross-sectional deformation on the properties of the en- 
capsulated molecules have also been suggested B121I13I1 . 

Among many kinds of intercalates, special attention has to be 
paid for dioxygen molecules (O2) confined into single-walled 
carbon nanotubes (SWNTs). An O2 molecule is the best known 
realization of a p-electron magnetic material; it has two un- 
paired electrons in a doubly degenerate antibonding orbital, 
thus being a spin S — 1 magnet in the ground state. In the 
bulk solid state, O2 assembly often exhibits anti-ferromagnetic 
ordering with mutually parallel arrangement jl4ll5ll . But con- 
finement into thin long cavity can lead ferromagnetic couplings 
of adjacent O2 molecules, in which the O-O bond orientations 
of two entities are orthogonal to each other In fact, 

the ferromagnetic ordering of a one-dimensional O2 chain is 
available through confinement into an armchair SWNT of the 
(6,6 ) chilarity, as evidenced by molecular dynamics simulations 
[18]. Structural control of such magnetic molecular aggregates 
is attractive in view of low-dimensional physicochemistry, and 
thus attempts based on various host materials have been made 
to stabilize the ferromagnetic order of absorbed O2 molecules 

A possible consequence of the ferromagnetic O2 molecular 
chain inside a SWNT is that the resulting magnetic field affects 



electronic properties of the host SWNT. Magnetic field appli- 
cation in general breaks time-reversal symmetry of the system. 
This symmetry breaking can lead feasible field-induced phe- 
nomena such as the Aharonov-Bohm effect |21, 23l and the 
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Landau level splitting I23I1 depending on the field strength and 
direction. Particularly in the 02-filled (6,6) SWNT, encapsu- 
lated O2 molecules are extremely close to the surrounding car- 
bon layer; therefore, their magnetic moments are thought to 
yield significant alteration in the quantum states of n electrons, 
which has remained unsettled so far. 

In this work, we theoretically investigate the electronic struc- 
tures of finite-length armchair (6,6) SWNTs encapsulating the 
ferromagnetic O2 chain. Spectral shifts caused by the O2 ab- 
sorption are found to be strongly dependent both on the tube 
length and the energy region to be studied. A signature of the 
O2 chain confinement manifests in the field-induced gap at the 
Fermi energy. Interestingly, the gap arises only when the tube 
length equals to a multiple of three in units of carbon hexagon, 
while it vanishes other tube lengths. This length-dependent gap 
formation may be observed in the electronic response of the 
02-filled nanotubes to alternate voltage with tens gigaheltz fre- 
quency. 

2. Model and Methods 

2.1. Hamitonian 

We consider the electronic structure of 02-filled (6,6) 
SWNTs utilizing the nearest-neighbor tight-binding model. 
The host graphitic sheet is the one rolled up in the direction of 
a C-C bond with the bond length of ac e = 0.142 nm, having 
a structural periodicity of the lattice constant of a — V3ac-c = 
0.246 nm along the tube axis. The tube length L is thus writ- 
ten by N y a, where N v is the number of hexagons aligned in 
the axial direction. Accordingly, 12N Y hexagons composed of 
24(N y + !)[= Nc] carbon atoms are involved in the system. The 
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Figure 1: Schematic of the Oi -filled armchair (6,6) SWNT with the tube ra- 
dius p = 0.407 nm. Small (colored in gray) and large (blue) circles are the 
constituent atoms of carbon and oxygen, respectively. Thick arrows indicate 
the magnetic dipole moments fio, associated with the unpaired electrons in O2 
molecules, which align with a separation rfo,-o 2 = 0-33 nm. The magnetic 
field profile B generated by single O2 molecule is depicted by curved paths. 

dangling bonds at the tube edges are assumed to be passivated 
by hydrogen atoms. 

The Hamiltonian matrix is built from the subspace spanned 
by the Nc wave functions of 2p z orbitals. The nearest-neighbor 
Hamiltonian is given by ( H — Yji,jJi,f\c j, where y h j is the 
transfer energy between the nearest-neighbor sites, and c.'(c,) 
is the creation (annihilation) operator of the n electron at the 
site i. The electronic states are obtained by diagonalizing the 
Nc x Nc Hermitian matrix. When the system is subjected to a 
magnetic field, the transfer integral y, j becomes y,-j = yoe 2m ^'' 
with yo = 2.7 eV, where faj = (e/h) f J A ■ dl is the Peierls 
phase and A is the vector potential corresponding to the mag- 
netic field. In the present system, the ferromagnetic O2 chain 
produces a spatially modulated magnetic field, and the mag- 
netic flux passing through each carbon hexagonal ring provides 
a finite value of 0y at edges of the hexagon. 

2.2. Magnetic field distribution 

The stable magnetic spin (S = 1) in an isolated O2 molecule 
generates the magnetic field Bo, at a point r = pe p + ze z in 
terms of the cylindrical polar coordinates (p, 0, z). The induced 
magnetic field is written by 

n jU 3r{r-n 0l )-r 2 no 2 

B 0l (r) = ? , (1) 

where fio 2 = -g e fi B S/fo with g e - 2, S = foe-, and p. B = 
efo/(2m e ); e and m e are the bare electron's charge and mass, re- 
spectively, and //o is the permeability of vacuum. Our task is to 
evaluate the normal component of Bo,, given by B ± = Bo, ■ e p , 
which exerts just on the graphetic sheet. Using the relations of 
f*o, -e p = and r-e p with p being the tube radius, the sum of the 
normal components, B^'(z), associated with all O2 molecules 



Figure 2: Spatial distribution of the normal component of the magnetic 
field caused by the O2 molecular chain in the hollow cavity of (6,6) SWNTs. 

encapsulated in the (6,6) tube is given by 

No, 

u _ -3f2 fi B y p(z - p - ndp 2 -o 2 ) 

2 * h\P^(z- P -ndo 2 -o 2 n 5ir 

Here, do 2 -o= 0.33 nm is the separation between adjacent O2 
molecules [18], and A^o, is the total number of the O2 molecules 
that position at z e [p,L - p] with equiseparation. We set p = 
0.407 nm according to the formula p - a V« 2 +m 2 + nm/ (2n) 
for nano tubes with (n, m) chirality. As a result, the magnetic 
flux O(z) that penetrates a hexagonal plaquette at z reads as 
<D(z) = V3Bf 1 a 2 /2, and the Peierls phase <pij can be deter- 
mined via the relation 2 hexagon <Pi,j = ^(z)/^0 with Oo = h/e 
being flux quantum. 

Figure [2] illustrates the spatial distribution of Bf'(z) for sev- 
eral tube lengths expressed by N y . Solid circles indicate the 
position of O2 molecules that may be contained in the host nan- 
otube. In the N y - 10 case, for instance, the six O2 to the most 
left (i.e., up to the O2 at z = 2.06 nm) can be included within 
the allowed region of z e [p,L - p] so that No, = 6. Simi- 
larly, No, = 10, 13 for /Vy = 15,20, respectively. We observe 
that B^'(z) oscillates with a slight amplitude in-between region 
but exhibits a large amplitude on the order of 20 mT near the 
tube edges. The amplitude suppression at the intermediate re- 
gion is attributed to cancelling out the oppositely-oriented field 
components associated with O2 molecules separated by ~ 1 nm 
each other. Neighborhoods of the tube edges are free from the 
suppression due to the absence of O2 both at z < p and z > L-p, 
and thus significant amount of the Peierls phase is imposed lo- 
cally to the electron's eigenmodes. As the scenario holds re- 
gardless of the tube length, the electronic states for longer (6,6) 
SWNTs are modulated to a degree by the O2 absorption, despite 
of slight magnetic moments of individual O2 molecules. 

3. Results 

Figure [3] displays the 02-induced shift in the electronic en- 
ergy spectra of short (6,6) SWNTs. The tube length is set to 



2 




• Ny=24 
■ Ny=25 

• Ny=26 



-1.5 -1.0 -0.5 0.0 0.5 1.0 
Pristine eigenfrequency E (THz) 



JJx10 



15 



10 



N 

X 

CD 



uj 5 



N 4 

a. 

LU -4 



Eg 



• Ny=24 
■ 25 

A 26 



-150 -100 -50 50 100 150 
E (THz) 



50 



100 

Ny 



150 



200 



Figure 3: Top: Spectral density of states of an empty (i.e., 09-unfilled) (6,6) 
nanotube. Bottom: O^-induced spectral shift for several tube lengths measured 
byiVy. 



Figure 4: HOMO-LUMO energy gap E„ as a function of the tube length Ny, E g 
grows monotonically with increasing N y and converges to ~ 13.6 GHz in the 
large N y limit. Inset: Enlarged view of the 02-induced eigenenergy shift in the 
vicinity of E = 0. 



be N y = 24, 25, 26 as simple examples, while much longer nan- 
otubes yield similar results as those in this plot. The horizon- 
tal axis E represents the eigenenergies of the pristine system 
(i.e., prior to O2 absorption), and the vertical axis AE gives 
their shifts caused by O2 absorption. Each discrete point in 
the plot indicates the quantum states of isolated finite-length 
(6,6) SWNTs, which are written by standing waves of the 
form i/f n (z) °c sin(k„z) with a discrete set of wavenumbers 
k„ = nn/(Nyd) where n = 1,2, ■• • ,N y Q. If the 2 chain 
were absent, AE = so that Nc discrete points would be line 
up horizontally in the bottom panel of Fig. [3] But in the current 
system, the C>2-induced magnetic field modulates the n elec- 
tron states, which results in upward or downward shift in the 
eigenenergies as quantified by AE in the figure. 

The data in Fig. [3] show pronounced shifts from the pristine 
values at discrete energies that correspond to the van Hove sin- 
gularities of the spectral density of states (top panel). A marked 
exception is the anomalous point-like shift at the band center 
(E = Ef — 0), which arises only in the case of N y = 24. 
This point-like energy shift implies an energy-gap formation 
between the highest occupied state and the lowest unoccupied 
state in the vicinity of E = 0. It should be emphasized that 
the energy gap no longer disappears in the other two cases of 
N y — 25, 26 as clarified visually in Inset of Fig. [4] 

We have revealed that the Cb-induced energy gap mentioned 
above is universally observed every when mod(N y , 3) = 0. Fig- 
ure|4]presents the N y dependence of the gap magnitude E g ; Note 
that only E g at N y = 3,6,9, •• • are plotted, because E g = 
whenever mod(N y , 3) ^ 0. The result exhibits a monotonic in- 
crease with N y followed by a convergence in the large N v limit 
with the asymptotic value of E g ~ 13.6 GHz. We mention that 
N y = 201 corresponds to the tube length ~ 60 nm, which scale 
is in a realm of the current fabrication techniques. It also fol- 
lows from Fig. [4] that the Zeeman splitting is irrelevant to the 
gap formation, since the energy scale fi^B in the present condi- 
tion is on the order of 0. 1 GHz, far less than the gap magnitude. 



4. Discussion 

The length-dependent gap formation at the Fermi level is un- 
derstood by considering the eigenlevel configuration for differ- 
ent N v . The left panel in Fig. [5] shows the energy dispersion 
relation for the infinitely-long (6,6) nanotube free from the fer- 
romagnetic O2 chain. There are fourteen branches reflecting 
the quantum confinement effect on the n electrons along the 
circumferential direction. Of these branches, ten thick curves 
are doubly-degenerate at every k points, while four thin curves 
are non-degenerate except for the crossing point at k — kp — 
2n/(3a). All the degeneracy is lifted when the ferromagnetic 
O2 chain is built inside the SWNT; it then leads the spectral 
shifts of all initially-degenerate eigenmodes that lie along the 
ten thick branches, which have been already demonstrated in 
the bottom panel of Fig. [3] 

The abovementioned scenario holds for finite-length (6,6) 
nanotubes, while the discreteness of allowed wavenumbers in 
the axial direction requires subtle modification. The right panel 
of Fig.[5]gives a drawing of the level configuration for A^, -long 
02-free (6,6) nanotubes in the vicinity of E = 0. Solid cir- 
cles indicate the eigenlevels for N y such that mod(A^ v , 3) = 0, 
and square ones for mod(A^ v , 3) + 0. In the former case, only 
one doubly-degenerate state (indicated by a large solid circle) 
resides just on the crossing point. In contrast, the latter case 
involves no degenerate state. Upon the O2 chain encapsulation, 
therefore, the degeneracy removal arises only in the A^,-long 
nanotubes with mod(N Y , 3) = 0. This mechanism accounts for 
the length-sensitive gap formation at the Fermi level. 

The degeneracy at the crossing point owes the availability 
of an appropriate value of the discrete wavenumber k„ = nAk 
with interval Ak = 7r/(N y d). When N y - 31 with integer £, the 
interval reads Ak = n/(3€a) and thus the discrete set of k n can 
take the value of the Fermi wavenumber kr = 2n/ (3a) at n = 2(. 
But when N y + 3(, the resulting interval Ak' = n/[(3{ ± l)a] 
forbids any k„ to coincides with kf since there is no integer n 
such that n/(3£± 1) = 2/3. 
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A^ v -long (6,6) SWNTs with N y being multiples of three; the gap 
vanishes for other N y . The length-sensitive gap formation is at- 
tributed to the removal of eigenlevel degeneracy at the Fermi 
level that arises only when mod(A^ v , 3) = 0. The gap magni- 
tude approaches 13.6 GHz in the long-A^ v limit, which suggests 
the possibility that the 02-induced gap can be confirmed in the 
alternative-current measurement of the host SWNT 
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Figure 5: Left: Energy dispersion relation for armchair (6,6) nanotubes free 
from the ferromagnetic O2 chain. Ten thick curves are doubly-degenerate at 
every k points, while four thin curves are non-degenerate except for the cross- 
ing point ax k = kf = 2w/(3o). Right: Eigenlevel configuration of an empty 
nanotube near the Fermi energy Ep = 0. Only when mod(A',,, 3) = 0, a doubly- 
degenerate state (indicated by a large solid circle) resides just on the crossing 
point. 



We suggest that the 02-induced gap with the magnitude of 
13.6 GHz can be detected through the dynamical conductivity 
measurement at low temperature. According to the Kubo for- 
mula 1 25, 2^], the dynamical conductance along the tube axis 



is found as 



G(co) = — V \J Z 



fu-fy 



fwj + Ejj — E v + id ' 



(3) 



where and f v are Fermi distribution functions and 6 is pos- 
itive infinitesimal. The magnitude of GicS) is primarily deter- 
mined by J~ v = (E M \J Z \E V ), in which J z is the z component 
of the current operator J = (-ie/h)[r, 'H] with a commutator 
[•"•, •*■]■ Suppose that ac electric voltage with frequency lo is 
applied to the 02-filled nanotubes. A sufficiently low temper- 
ature of k^T < Hot, Eq. © is reduced to G{o>) cc \J + _\ 2 /u>, 
where /+_ = (E+\J\E-) and £+(-) is the eigenlevel immedi- 
ately above (below) Ep = 0. In the present system, /+_ has 
a modestly large value, which follows from the numerical re- 
sult of |2/I C i {E+\ri){ri +e z \E-)\ ~ 0.2 under the normalization 

condition of Yjf~i \i r i\E+(-))\ 2 = 1. Hence, the dynamical con- 
ductivity will exhibit a sharp peak at u> ~ 13.6 GHz in a low- 
temperature measurement at less than 1 K. A large-scale com- 
putation technique for the Kubo-formula II27[ 12811 incorpolated 
with ab initio calculations 1E9I1 enable a quantitative examina- 
tion of the conjecture, which we will perform in the future work. 

5. Summary 

We have theoretically considered the electronic spectral shift 
of the (6,6) SWNT encapsulating a linear ferromagnetic O2 
chain into the hollow cavity. O2 absorption causes an energy 
gap between the highest and lowest occupied states only in the 
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